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ABSTRACT 

The  mechanical  reliability  of  long  lengths  of  glass  fiber  optical 
waveguides  is  important  to  the  development  of  optical  communications 
systems.  In  the  present  work,  long  silica  fibers  were  drawn  and  coated 
in-line  with  polyethylene.  The  strength  of  specimens  from  0.05  to 
11.96  meters  was  measured,  and  the  strength  distribution  of  119 
specimens,  1.06  meters  long  was  also  measured. 

Three  methods  of  estimating  the  parameters  of  the  Weibull 
statistical  distribution  are  compared,  and  the  influence  of  bimodal 
populations  on  the  shape  of  an  assumed  Weibull  distribution  are 
explored.  A technique  for  extracting  a bimodal  distribution  is  demon- 
strated and  applied  to  the  test  sample  of  119  silica  fibers. 
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1 INTRODUCTION 

Optical  waveguides  using  long  glass  fibers  about  100  to  200  pm 
diameter  require  reliable  and  predictable  strength  characteristics 

I in  service.  The  problem  is  basically  that  of  predicting  the  probab- 

ility of  failure,  or  the  minimum  extreme  strength,  of  a one  kilometer 

: glass  fiber  from  mechanical  strength  measurements  on  much  shorter 

samples.  The  nature  of  brittle  failure,  the  need  to  include  a 

I 

significant  size  parameter  in  scaling  up  the  test  by  a factor  of  1000 
or  more,  and  the  need  to  predict  the  population  behavior  from  a limited 
sample  all  require  the  application  of  statistical  methods  of  design 
reliability. 

The  purposes  of  the  research  carried  out  under  this  task  were  the 
following:  to  develop  techniques  for  drawing  long  silica  fibers  with 

; 

I 

i in-line  pol)nner  coating, to  test  the  mechanical  strength  of  long  glass 

i 

I fibers, and  to  examine  the  validity  of  techniques  commonly  employed  to 

i 

I apply  the  Weibull  statistical  model  to  brittle  strength  data. 

[ The  experimental  result  of  this  work  has  already  been  summarized 

in  the  Final  Report,  "Optical  Coupling  Techniques,  April  1,  1976.(1) 
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II  FIBER  FABRICATION  AND  STRENGTH 


1.  Drawing  Silica  Fibers 

The  final  fiber  drawing  device  developed  for  the  production  of 
coated  silica  fibers  is  shown  schematically  in  Fig.  1.  From  the  top 
downwards,  the  following  features  were  incorporated  in  the  apparatus: 

(a)  Variable  speed  feed  mechanism  to  lower  the  preform  into  the  furnace 
at  a constant  rate.  (b)  Graphite  sucecptor  nitrogen  protected  induction 
furnace  capable  of  2000‘‘C,  shown  in  Fig.  2.  Milmaster  Model  SSE-5R  to 
monitor  a preset  diameter  between  100  to  200  pm  within  5%.  (d)  Polyethy- 

lene coating  unit  consisting  of  a pressurized,  heated  annular  tank,  shown 
in  Fig.  3.  (e)  A cooling  station  for  air  cooling  or  water  mist  cooling 

of  the  polyethylene  coating.  (f)  Precision  machined  aluminum  winding 
drum  with  variable  speed  drive,  (g)  Traversing  base  driven  from  the 
winding  motor  to  maintain  constant  pulling  position.  This  is  the  same 
apparatus  used  in  pulling  special  waveguide  shapes  for  optical  coupling 
structures.  (2) 

The  silica  fibers  used  in  this  study  were  drawn  from  as  received 
5 mm  diameter  silica  rods.*  They  were  not  duplex  core  cladding  fibers 
nor  were  optical  transmission  characteristics  controlled  or  measured. 
Typical  diameters  were  125  to  140  ym  with  a coating  thickness  of  125  ym. 

2.  Strength  Tests 

Attempts  to  test  the  strength  of  fibers  several  hundred  meters 
long  used  a drum  to  drum  device  described  earlier.  (2)  This  approach 
was  abandoned  because  of  occasional  random  failure  of  the  fiber  while 


J 

1 

I 


k 


resting  under  tension  on  the  winding  drum  and  the  uncertainty  and 
*0ptosll  Til  from  Amersll,  Inc.,  Hillside,  N.J. 


VARIABLE  SPEED  FEED 


SILICA  ROD  PREFORM 


INDUCTION  FURNACE  (FIG.  2) 


DIAMETER  MONITORING 


POLYETHYLENE  COATING  (FIG.  3) 


•COOLING  STATION 


VARIABLE  SPEED 
WINDING  DRUM 


■TRAVERSING  BASE 


Fig.  1.  Fiber  Drawing  Apparatus.  Basic  design  after  Bell  Laboratories 
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non-uniformity  true  stress  introduced  by  curvature  on  the  drums.  The 

best  estimate  of  minimum  strength  found  in  fiber  lengths  up  to  300  m 
2 

was  66. A MN/m  (9500  psi) . 

Fibers  8 to  12  meters  long  were  tested  while  suspended  vertically 
in  a utility  well.  Constant  loading  rate  was  obtained  by  running 
water  into  a container  from  a constant  head  supply.  Individual  results 
from  these  tests  are  given  in  Table  I. 

Fibers  l.A  meter  and  less  in  length  were  tested  in  an  Instron 
Testing  machine  at  constant  strain  rate  which  is  equivalent  to  constant 
load  rate  for  materials  in  their  elastic  range.  Table  II  gives  strength 
results  for  various  size  fibers  while  Table  III  gives  results  for  a 
large  sample  of  119  fibers,  1.06  meters  long. 

A major  experimental  difficulty  throughout  this  program  was  that 
of  obtaining  effective  mechanical  gripping  through  or  on  the  polyethylene 
coating.  The  results  reported  on  Tables  I,  II  and  III  used  grip  pads 
made  from  polyethylene  blocks  which  were  fused  to  the  ends  of  the  fiber 
with  molten  polyethylene.  This  was  a very  time  consuming  and  tedious 
procedure.  At  the  high  stress  levels  commonly  obtained  for  short  gage 
length  fibers  (Table  II)  many  failures  occurred  by  shearing  of  the 
coating,  and  fiber  strengths  greater  than  about  2750  MN/m2  (AOO  ksi) 
were  very  difficult  to  measure.  Some  later  trials  using  roller  grips 
where  the  fiber  was  wrapped  several  turns  around  a 3 cm  diameter  roller 
surfaced  with  rubber  pads  appeared  to  be  a promising  technique. 

3.  Sequential  Tests  and  the  Weakest  Link  Theory 

A basic  postulate  in  the  theory  of  brittle  fracture  is  that  failure 


occurs  at  the  largest  flaw  in  the  system  - the  so-called  Weakest  Link 


Table  I Tensile  Strength  of  Coated  Silica  Fibers. 

Test  No.  Gage  Length  Fracture  Strength 

2 

meters  MN/m 


2 

11.96 

551.6 

3 

11.96 

537.8 

4 

11.99 

448.2 

5 

8.38 

413.7 

6 

8.28 

524.0 

7 

8.33 

475.8 

9 

8.31 

620.5 

10 

8.41 

572.3 
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Table  II  Tensile  Strength  of  Coated  Silica  Fibers 
Tested  in  an  Instron  Testing  Machine 


Test  No. 

Gage  Length 
meters 

Fractun 

MN/m"^ 

1 

0.05 

>3102.7 

2 

It 

II 

3 

tl 

>2413. 

4 

ft 

>2068 

5 

tl 

>2758 

6 

ft 

>3447 

7-17 

0.076 

>1930 

8-18 

II 

>2102 

9-22 

It 

>1220 

10-23 

tl 

>1772 

11-24 

It 

882.6 

12 

0.102 

>2758 

13 

0.254 

>3102 

14 

It 

>2758 

15 

tl 

>3102 

16 

II 

>3448 

17 

It 

>3448 

18 

If 

>3448 

19 

tl 

>2758 

20  . 

II 

>2758 

21 

0.333 

>3102 

22-15 

0.356 

>1448 

22-16 

II 

1930.6 

23 

1.09 

1551.4 

24 

1.14 

482.6 

25 

II 

>1379 

26 

tl 

2413.2 

27 

1.17 

620.5 

28 

1.19 

827.4 

29 

1.22 

482.6 

30 

It 

>2413 

31 

1.27 

413.7 

32-20 

1.37 

958.4 

33-21 

II 

827.4 

Strengths  with  > symbol  indicate  failure  shear  of  the 


Table  III.  Ordered  Strength  Results  for  119  coated  silica  fiber 
Specimens  1.06  m long. 


Class  Strength  Frequency 

Number  MN/m^ 


Cumulative  Probability 

Frequency  (1/n+l) 


1 

24m.? 

1 

1 

.000 

? 

337. Q 

1 

2 

.017 

4HP.4 

1 

3 

.025 

u 

SIO.? 

1 

4 

.033 

s 

530, Q - 

1 

5 

.042 

537, ft 

1 

6 

.050 

7 

5<3«).Q  - 

1 

7 

.058 

H 

627.4 

1 

0 

.067 

Q 

634.3  

2 - 

10 

.0«3 

m 

64ft. 1 

1 

11 

.092 

n - 

661. Q 

1 

. 12 

.100 

1? 

696,4 

1 

13 

.108 

1 

710.?  - - 

1 

14 

.117 

M* 

737, ft 

1 

15 

.125 

IS 

744.7-- 

1 

16 

.133 

I#' 

75ft. S 

1 

17 

.142 

1 -f 

77Q  _ 1 

1 

- 18 

.150 

786.0 

2 

20 

.167 

IQ 

79?. 9 _ 

? 

- ?2 

.183 

yn 

799. ft 

2 

?4 

.200 

P 1 

MnK^7 

. P - — . 

. - . ?6  - 

.217 

?? 

020. S 

1 

27 

.225 

PP7  'i 

1 

' ?ft 

.233 

?u 

034,3 

2 

30 

.250 

PC 

AA 1 P 

- . . - 34 

.263 

840.1 

1 

35 

.292 

P7 

B5S.  0 

. 1 - 

36 

.300 

060,0 

i 

37 

.308 

PQ 

88?. f,  

2 

. . . 39 

.325 

•»n 

903.2 

1 

40 

.333 

31 

910.1-  -- 

2 

-42 

.350 

923.9 

? 

44 

.36  7 

n 

937.7 

. ? 

j46 

.383 

34 

944.6 

1 

47 

.39? 

3S 

<15 1 .5 

^ . 

AQ  _ 

.408 

If. 

958.4 

i 

50 

.417 

37 

96S  ■ ■* 

1 

*;l 

.425 

3h 

97?.? 

i 

52 

.433 

3Q 

979.1 

1 

53 

.442 

40 

986.0 

1 

54 

.450 

99?. Q 

.4 

- 58  - 

.463 

to? 

999. ft 

1 

59 

.492 

43 

. 1006.7 

1 . 

. 60 

.500 

44 

1013.6 

2 

62 

.517 

toS 

1027.4 

1 

63 

.525 

9A 


Tabic  III.,  (continued) 


Class 

Strength 

Frequency 

Cumulative 

Probability 

Number 

MN/m^ 

Frequency 

(i/n+1) 

4A 

1034.? 

2 

65 

.542 

U1 

1043.0  

1 

66 

.550 

UA 

1054.4 

2 

68 

.567 

44 

107S.f. 

2 

70 

.583 

SO 

1089.4 

1 

71 

.592 

SI 

1110.1 

1 

72 

.600 

5? 

1117.0 

1 

73 

.608 

s:?  - 

1123.4 

2 

75 

.625 

S4 

1 144.3 

2 

77 

.642 

ss 

1 151. S 

1 

78 

.650 

Sf* 

1 158.4 

1 

79 

.658 

S7 

1172.? 

J 

80 

.667 

SB 

1179.0 

2 

82 

.683 

SQ 

119?. 8 

2 

. 84 

. 700 

hO 

1 199.7 

1 

85 

.708 

SI 

12^1.1 

1 

86 

.717 

S? 

1254.4 

« • 

87 

.725 

hi 

1222.5 

88 

.733 

1303.? 

1 

89 

.74? 

hS 

1316.9 

1 - « ... 

- 90 

.750 

hh 

1330.7 

1 

91 

.758 

h7  - 

1365.? 

. .....  J 

92 

.767 

hft 

14?0.4 

1 

43 

.775 

hO 

1441.1  

- - 1 . - — 

94 

.783 

70 

1 454 . 8 

1 

45 

.792 

. 74 

1475. S 

- - 2 

47 

anQ 

7? 

148?. 4 

2 

99 

.825 

71 

- 1496.?  — 

1 

- 100 

.833 

74 

1537.6 

1 

lOl 

.842 

7S  _ 

— - 167S.5 

1 

102 

.850 

7h 

17?3.7 

i 

103 

.058 

77 

1751,3-  

1 n/. 

.867 

7fi 

1806.5 

1 

105 

.875 

_ 74  -- 

1847.4 

4 - 

- 106 

.883 

flO 

1972.0 

1 

107 

.892 

B4  . 

- 1992.7 

J 08 

.400 

fi? 

2144.3 

1 

109 

.908 

- a.i 

?324,h 

110 

.417 

84 

?427.0 

1 

. Ill 

.925 

8S- 

2463. 4 

f 

112 

.933 

8h 

2640.8 

1 

113 

.442 

87 

V716.6 

. - . — 1 - - — . 

114 

.950 

88 

2875.2 

115 

.958 

84 

- 2978.6  . 

116 

.967 

40 

3385.4 

1 

117 

.975 

91 

358S.4  

• I . 

- 118 

.983 

4? 

3771 .6 

1 

119 

.992 

.1 — t, 
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Theory.  In  order  to  test  this  hypothesis,  broken  sections  of  longer 
fibers  were  retested.  The  results  are  shown  in  Table  IV.  In  every 
case,  the  secondary  fractures  were  at  higher  stresses  than  the  initial 
fracture. 

Not  only  does  this  agree  with  the  weakest  link  theory,  but  also 
shows  that  the  polyethylene  coating  was  effective  in  preventing 
externally  caused  damage  to  the  fibers  during  handling  and  testing, 
and  that  mechanical  properties  of  the  fibers  were  not  degraded  by 
subjecting  them  to  prior  stressing.  These  latter  points  are  important 
considerations  in  proof  testing  for  fiber  reliability  and  degradation 
of  fiber  strength  in  service. 

A.  Statistical  Analysis  of  Strength  Results 

The  data  of  Table  III,  119  specimens,  1.06  meters  long,  was 
analyzed  in  detail  to  determine  if  the  two-parameter  Weibull  distribu- 
tion would  adequately  model  these  results  and  permit  valid  inferences 
on  the  effect  of  specimen  size  on  strength.  In  the  following  sections, 
the  Weibull  model  and  methods  of  estimating  Weibull  parameters  are 
briefly  reviewed.  The  effects  of  variations  of  Weibull  parameters 
on  the  shape  of  calculated  distribution  functions  are  shown,  and  the 
Influence  of  blmodal  distributions  on  the  commonly  used  linear  plotting 
technique  are  calculated.  The  data  of  Table  III  is  then  re-considered 
as  a bimodal  distribution,  and  the  size  effect  on  strength  is  plotted. 


label  IV  Sequential  Strength  Tests  of  Silica  Fibers 


FIRST  TE 

ST 

SECOND  TEST 

THIRD  TEST 

Gage 

Strength 

Gage 

Strength 

Gage 

Strength 

Length 

MN/m^ 

Length 

MN/m2 

Length 

MN/m^ 

meters 

meters 

meters 

Test  9,  Table  I 

8.33 

620.5  ~ 

882.6 

703.3 

675.7  

>-  0.04 

855 

0.08 

>1710 

0.08 

>2758 

0.30 

>1517 

0.30 

>1237 

Test  32-20,  Table  II 


1.37 

958.4  — 

0.06 

1.02 

1.22 

>2517 

>2034 

1172. 

Test  7, 

Table  I 

8.33 

475.8  

0.91 

827. 

1.07 

586.1 

1.07 

599.9 

1.07 

910.1 

1.07 

792.9 

1.07 

861.9 

1.52 

620.5 

Ill  STATISTICS  OF  THE  STRENGTH  OF  BRITTLE  MATERIALS 


The  mechanical  failure  of  a brittle  material  is  caused  by  the 
propagation  of  a microscopic  crack  or  Griffith  flaw.  In  the  case 
of  glass  fibers,  these  flaws  are  a few  micrometers  or  less  in  size 
and  occur  primarily  on  the  surface  of  the  material.  The  flaw  size, 
the  number  of  flaws  and  their  frequency  of  occurrence  may  be  de- 
scribed by  statistical  distributions.  The  general  problem  is  one 
of  predicting  the  minimum  extreme  value  of  strength  resulting  from 
a sample  (flaw)  population  of  a given  size.  The  exact  solution 
for  extreme  value  distributions  has  been  given  by  Gumbel  (3,4) 
and  Epstein  (5,6)  has  published  extreme  value  formulations  for  a 
number  of  different  statistical  distributions.  One  type  of  extreme 
value  model  which  has  wide  applicability  in  failure  statistics  of 
many  kinds  and  brittle  failure  in  particular  is  the  Weibull  model  (7,8). 

The  Weibull  cumulative  distribution  function  (CDF)  relating 
the  probability  of  failure,  F(x)  to  the  strength,  x is 

F(x)  = 1 — exp  {-(x/o)^}  (1) 

where  a and  X are  the  scale  parameter  and  the  shape  parameter  re- 
spectively.* The  corresponding  probability  density  function  (PDF) 
defined  as  dF(x)/dx  is 

f(x)  * (X/o) (x/o)^“^  exp{-(x/o)^}  (2) 

0<x,a,X 

Transposing  and  taking  the  logarithm  twice  transforms  Eq.  1 into 

In  In  (1/(1-F(x)))  - Xln  x - Xlna  (3) 


*The  notation  and  nomenclature  will  follow  that  of  K.V.  Bury 

Statistical  Mode  Is  in  Applied  Science.  Ref.  9,  where  the  deriva- 
tion and  properties  of  this  and  other  models  are  given  In  detail. 
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. The  effect  of  specimen  size  on  observed  strength  is  actually 

observed  as  the  variation  of  the  extreme  value  of  flaw  severity 
when  the  flaw  sample  size  is  Increased,  l.e.  the  number  of  flaws 
"sampled"  in  a given  individual  strength  test.  Following  Bury  (9) 
the  sample  size  enters  the  Weibull  CDF  as 

F(x)  = exp  {-  (L/Lj^)  (x/o)^}  (4) 

where  is  a unit  length  (or  area  or  volume)  containing  flaws 
with  associated  content  a . L is  some  greater  length  with  the 
same  number  of  flaws  per  unit  area  and  a=aj^/(Nj^) 

At  constant  probability  of  failure,  e.g.  0.5,  the  relation  of  size 
to  failure  stress  is 

(5) 

This  was  originally  established  by  Weibull  (7)  on  empirical  grounds 
and  utilized  recently  by  Tarliyal  and  Kalish  (10),  Kalish  et  al.(ll), 
^ ' Kurkjian  et  al.  (12),  and  Mauer  et  al.  (13). 
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IV  ESTIMATION  OF  THE  WEIBULL  PARAMETERS  a and  X 
1.  Linear  Plotting 

The  most  common  technique  for  estimating  the  parameters  o and 
X is  to  plot  the  observed  data  in  the  linear  form  of  Eq.  3.  The 
expected  value  of  F(x)  is  obtained  from  the  ordered  data  where 
F(x)  “ i/(n+l)  (6) 

in  which  i is  the  1th  order  of  failure  and  n is  the  total  number 
of  specimens  in  the  sample.  In  addition  to  initial  estimates  of 
o and  X,  this  plot  also  gives  information  as  to  the  suitability 
of  the  proposed  model.  Substantial  systematic  curvature  in  what 
should  be  a linear  plot  indicated  the  need  for  a third  parameter 
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In  the  Weibull  function  (9)  or,  as  will  be  showm  later,  a blmodal 
or  multi-modal  population. 

2.  Maximum  Likelihood  Method 

A second  method  called  Maximum  Likelihood  (ML)  minimizes  Che 
variance  of  the  parameters  a and  X,  Thoman,  Bain  and  Antle  (14) 
analyzed  this  technique  for  the  Weibull  distribution  and  published 
confidence  Intervals  for  the  parameters  and  unblaslng  factors  for 
the  shape  parameter  as  a function  of  sample  (test)  size.  The 
Maximum  Likelihood  value  of  A is  found  by  an  interative  procedure 
and  the  expression  given  by  Thoman  (Ref.  14,  Sec.  5)  converges  to 
within  + 0.001  in  about  four  iterations  using  starting  values 
from  the  linear  plot  of  Eq.  3.  The  ML  estimate  of  o is  then  found 
directly. 

3.  Non-Linear  Least  Squares  Estimation 
Davies  (15)  proposed  using  direct  non-linear  least  squares 
curve  fit  to  the  CDF,  Eq.  1 as  an  improvement  in  some  situations 
over  linear  fits  to  Eq.  3 or  maximum  likelihood  estimates.  Non- 
linear curve  fitting  programs  are  not  generally  available,  how- 
ever, a well  documented  public  computer  program  called  NLWOOD* 
has  been  described  by  Daniel  and  Wood  (16) . It  Is  an  Iterative 
routine  which  will  evaluate  a and  X directly  from  Eq.  1.  While 
the  program  does  not  minimize  the  squared  residuals,  it  permits  a 
user  selected  criteria  of  relative  change  in  fitted  parameters 
and/or  sum  of  squared  residuals  as  a test  of  fit.  It  provides 
confidence  intervals  on  the  parameters,  sum  of  the  squared  resi- 
duals and  useful  graphical  Information  on  the  individual  residuals. 

A comparison  of  the  three  methods  for  estimation  of  Weibull  para- 
meters has  also  been  given  by  Heavens  and  Murgatroyd  (17). 


*The  source  program,  FORTRAN  listing  and  User's  Manual  is  available 
from  the  SHARE  Library  (Number  360D-13. 6. 007) , for  IBM,  and  the 
VIM  Library  (Number  G2-CAL-NLW00D) , for  CDC. 


V WEIBULL  PARAMETERS  AND  THE  RESULTING  DISTRIBUTING  FUNCTIONS 


Theoretical  Weibull  distributions  were  generated  by  assuming 
values  for  o and  X,  and  a sample  size  of  60  and  calculating  the 
corresponding  strength  using  Equations  1 and  6.  Plots  were  made 
of  Equations  1,  2 and  3 for  the  assumed  parameters  to  observe  the 
effect  of  systematic  parameter  variations  on  the  shapes  and  loca- 
tions of  the  various  curves.  In  addition,  the  effect  of  bimodal 
distributions  on  the  shape  of  "linear"  plots  was  studied  by  merg- 
ing two  sets  of  calculated  strength  values,  reordering,  and  re- 
plotting  the  merged  set  as  if  it  was  a single  unkown  sample. 

The  Weibull  distributions  selected  were  similar  to  distributions 
observed  in  our  laboratory  for  the  strength  of  7 and  10mm  glass 
rods  in  3-point  bending. 

Figure  4 shows  the  effect  of  varying  a with  constant  X, 
(Table  V).  This  is  essentially  the  size  effect  for  a given 
population  where  the  smallest  specimen  size  is  represented  by 
curve  D and  the  largest  size  by  curve  A,  nearly  2200  times  larger 
than  D.  (An  extrapolation  of  50cm  fiber  tests  to  1km  is  a size 
effect  of  2000).  Although  the  Weibull  parameter  X determines  the 
slope  of  the  linear  plot  (Fig.  4c)  it  is  clearly  incorrect  to  say 
that  X determines  the  "steepness"  of  the  CDF  (Fig.  4a)  or  the 
dispersion  of  the  PDF  (Fig.  4b).  These  two  features  which  relate 
to  the  variability  of  a set  of  measurements  around  the  mode  or 
average  are  a function  of  both  the  shape  parameter  and  scale  para- 
meters and  only  in  the  case  where  the  scale  parameter  (or  the  more 
directly  measured,  average  strength)  is  the  same  will  the  shape 
parameter  alone  determine  the  "steepness"  of  the  CDF. 


Figs.  4a,  4b.  Calculated  Welbull  CDF  and  PDF  curves  with  X constant 
and  o increasing  left  to  right.  See  Table  V.  This  is  equivalent 
to  change  in  specimen  size  (e.g.,  length)  in  a constant  population 
in  which  A specimens  are  nearly  2200  times  longer  than  E specimens. 
The  PDF  curves  are  scaled  to  equal  areas. 
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Fig.  Ac.  ‘ Linear  forms  of  the  Weibull  functions  of  Aa  and  Ab.  The 
data  points  were  calculated  from  two  individual  functions  and  then 
merged,  reordered  and  plotted  as  if  coming  from  a single  sample. 


TABLE  V 

PARAMETERS  USED  IN  CALCULATING  WEIBULL  FUNCTIONS 

Curve 


A 

B 

C 

D 

E 

,Flg.  A 

a 

200 

300 

AOO 

500 

600 

X 

7.0 

7.0 

7.0 

7.0 

7.0 

Fig.  5 

0 

200 

300 

AOO 

500 

600 

X 

5.0 

7.0 

9.0 

11.0 

13.0 

Fig.  6 
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AOO 
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AOO 

AOO 
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• The  points  in  Figure  4c  are  the  merged,  calculated  values 
from  curves  A-B,  A-C,  A-D,  and  A-E,  i.e.  sampling  values  which  might 
arise  if  the  real  populations  were  an  unknown  bimodal  distribution. 
As  seen  in  Figure  4c,  X for  lower  strength  distribution  (A)  would 
be  fairly  accurately  determined  from  the  slope  of  the  low  strength 
points.  However,  Ina,  which  is  the  value  on  the  abscissa  where 
the  linear  plot  intersects  the  horizontal  line  of  In  In  (<}))  = 0, 
would  have  substantial  error.  Virtually  no  accurate  information 
can  be  obtained  from  this  plot  for  the  higher  distribution  of  the 
pair. 

Figure  5 is  an  example  where  both  a and  X Increase  (Table  V) . 
This  produces  PDF  curves  with  very  similar  central  portions  but 
different  tails  skewed  to -the  left.  Again  errors  in  estimating 
0 from  the  merged  samples  will  arise  if  the  lower  strength  popula- 
tion is  assumed  to  coincide  with  the  data  points.  This  error 
occurs  in  the  lower  part  of  Figures  4c  and  5c  because  the  order 
plotting  function  for  each  of  the  true  populations  is  i/(n+l) 
while  for  the  merged  population  it  was  i/(2n+l).  As  pointed  out 
by  Tarlyal  and  Kallsh  (10), to  separate  bimodal  distributions,  one 
should  know  the  fraction  of  observations  belonging  to  each  dis- 
tribution. With  this  information  and  only  the  lower  part  of  the 
strength  data,  e.g.  30  of  the  60  values,  the  order  plotting  would 
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be  correct  and  more  accurate  estimates  of  o and  X would  result. 
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VI  ANALYSIS  OF  E}J>ERIMIiNTAl.  STRENGTH  DATA 

Figure  7 shows  the  direct  CDF  plot  for  the  119  specimens  1.06m 
long  listed  in  Table  III.  The  dashed  line  is  the  curve  calculated  using 
the  ML  estimates  of  the  Weibull  Parameters,  and  the  solid  line  was 
obtained  using  parameters  from  the  non-linear  direct  fit.  The  deviation 
at  the  higher  strength  is  very  similar  to  that  seen  by  Davies  (15)  in 
silicon  nitride.  Figure  8 shows  the  "linear"  plot  of  this  data  assuming 
Equations  3 and  6. 

One  possible  correction  for  the  curvature  in  Figure  8 is  the 
Introduction  of  a location  parameter  p in  the  Weibull  CDF  as 

F(x,u,o)  = 1 - exp  {-(^~)  (7) 

This  corresponds  to  a minimum  strength,  y,  below  which  failure  does 

2 

not  occvr*  A third  parameter  was  introduced  in  steps  of  35  MN/m 

2 

from  zero  to  the  minimum  observed  strength  of  248  MN/m  , and  the 
resulting  CDF  curve  calculated  at  each  trial.  Although  the  ML  CDF 
fit  (as  well  as  the  non-linear  direct  fit)  improves  slightly  as 
judged  by  a reduction  in  sum  of  squared  residuals,  the  curvature 
of  the  linear  plot  or  the  fit  at  the  upper  end  of  the  CDF  was  not 
significantly  Improved. 

Based  on  Fig.  6 and  other  similar  trials,  the  curvature  in 
Fig.  8 was  assumed  to  arise  from  an  overlapping  bimodal  distribu- 
tion. The  parameters  of  the  two  distributions  were  estimated  as 
follows:  Based  on  straight  line  portions  at  the  ends  of  the  curve 
of  Fig.  8,  the  population  totals  were  estimated  to  be  80  specimens 
for  the  lower  portion  and  39  for  the  upper  portion.  These  sample 
sets  will  be  called  Set  1 and  Set  II  respectively.  Using  these 
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population  sub-sets  and  assuming  the  lower  40  and  upper  30  data 
points  of  the  total  set  to  belong  to  each  of  the  Individual  dis- 
tributions alone,  the  30  data  points  were  re-ordered  against  their 
respective  sub-set  population.  Welbull  parameters  were  estimated 
using  either  straight  line  probability  plotting  or  non-linear 
curve  fitting  to  this  partial  data.  Maximum  likelihood  estimates 
could  not  be  used  because  the  30  points  are  not  the  entire  sample. 

Two  pairs  of  Welbull  parameters  o and  X were  estimated  from 
these  plots.  Theoretical  sets  of  80  and  39  data  points  were  then 
calculated  using  these  parameters,  and  two  theoretical  sets,  1 and 
11  were  merged  and  re-ordered  to  correspond  to  the  original  experi- 
mental observation.  The  ordered  positions  of  Individual  points 
belonging  to  set  1 and  Set  11  were  determined  In  the  merged 

population.  For  example,  in  the  merged  119  theoretical  points,  the 
following  positions  were  occupied  by  points  from  Set  11:  Position 
No.  l,4,8,14,21,30,39,49,59,67,75,82,87,91,93-to-119.  These  same 
positions  were  then  selected  from  the  ordered  experimental  data 
and  assigned  to  Set  II  with  the  remaining  assigned  to  Set  I. 

Welbull  parameters  were  then  recalculated  for  the  experimental  data 
and  the  results  are  given  in  Table  VI  and  Fig.  9.  There  is  a 
large  variation  in  the  parameters  estimated  by  the  three  different 
techniques  for  the  full  set  of  119  specimens.  The  estimation 
techniques  are  much  more  consistent  for  the  assumed  subsets. 

One  can  now  calculate  the  influence  of  these  two  distributions 
on  the  strength  of  long  fibers,  e.g.  1000  m.  Using  Eq.  4,  the 
failure  probability  for  each  distribution  at  1km  and  various  stresses 
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is  obtained.  The  total  failure  probability  is  then 
^’lotal  = 2 Pj^/3  + Pj.j/3 

assuming  the  fraction  of  the  populations  from  each  distribution 
remains  2/3  and  1/3.  The  results  are  given  in  Fig.  11  which  shows 
the  calculated  size  effect  on  strength  using  Eq.  1 and  the  joint 
probability  of  two  distributions.  Data  points  for  the  experimental 
results  for  fiber  lengths  from  1.06  to  11.93m  are  also  shown. 

In  calculating  the  contribution  to  failure  from  the  two  distribu- 
tions it  was  found  that  the  tall  of  the  low  A distribution  dominated 
at  long  lengths.  It  is  interesting  that  the  data  from  which  this 
distribution  was  extracted  was  in  the  higher  strength  portion  of  the 
original  experimental  measurement. 


Fig.  11.  Average  strength  of  experimental  fibers  and  the 
calculated  0.5  probability  from  the  joint  distribution. 


VII  CONCLUSIONS 


As  In  any  statistical  modeling,  the  foregoing  analysis  Is  an 
attempt  to  define  functions  which  represent  the  underlying  phenomenon 
The  fact  that  one  can  obtain  functions  which  appear  to  do  this  Is  not 
a confirmation  of  the  nature  of  the  underlying  phenomenon  but  with 
hypothesis  testing,  only  places  some  confidence  measurement  on  the 
validity  of  the  assumed  model. 

In  mechanical  property  testing  of  brittle  materials  where 
blmodal  distributions  are  suspected.  Independent  physical  confirma- 
tion of  different  failure  modes  and  measurement  of  the  fraction 
belonging  to  each  distribution  should  be  attempted.  Direct  fracto- 
graphlc  analysis  which  would  Indicate  different  types  of  flaw 
origins  would  be  useful.  However,  It  would  be  very  tedious  and 
perhaps  Impossible  for  large  numbers  of  high  strength  fibers.  The 
analytical  techniques  recently  published  by  Matthews  et-  al.  (19) 
may  be  useful  In  defining  different  Independent  populations  of  flaw 
origins . 

Maximum  Likelihood  method  of  estimation  of  Welbull  parameters 
is  the  most  attractive  since  no  ordering  of  the  data  Is  required, 

Che  convergence  Is  very  rapid,  and  the  form  of  the  expression 
given  by  Thoman  et  al.  (14)  easily  adapted  to  a computer.  Thoman 
et  al.  have  also  provided  the  analysis  of  the  Inferences  on  Che 
Welbull  parameters.  However  the  ML  method  can  not  give  correct 
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results  If  the  underlying  distribution  is  not  in  fact  a single 
Welbull  distribution.  For  this  reason,  probability  plotting  should 
be  employed  as  a check  on  the  assumption  of  the  statistical  model 
and  to  furnish  a first  estimate  of  \ needed  for  the  ML  iteration. 

When  the  estimates  obtained  by  least  squares,  non-linear  direct 

fit  and  Maximum  Likelihood  are  widely  different,  it  is  an  indication 

that  the  assumption  of  a single  modal  Welbull  distribution  is  Incorrect. 
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